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The major goal of the high-precision studies of ro- vibrational states in the hydrogen molecular ions 
is to provide an alternative way for improving the electron-to-proton mass ratio, or the atomic mass 
of electron. By now the complete set of relativistic and radiative corrections have been obtained for 
a wide range of ro- vibrational states of and IfD"'' up to order RooCt'^- In this work we complete 
calculations of various contributions to the RocOt^ order by computing the relativistic corrections to 
the binding energy of electron. 

In recent years several laser spectroscopy experiments have been proposed [1, 2] for high precision measurements of 
the vibrational spectrum of the hydrogen molecular ions Hj" and HD"*". These experiments have metrological interest 
and are aimed at a sup-ppb precision. In order to improve the present accuracy of the electron-to-proton mass ratio 
[3] the uncertainty of the spectroscopic data (as well as of the theoretical calculations of the spectra to compare to) 
should be below 1 part per billion(l ppb). To meet these stringent requirements, the theoretical calculations should 
achieve at least a level of 10 kHz (or ~10~^^ in atomic units). 

While the variational calculations of the nonrelativistic energies have reached a numerical precision of 10"^*^ — 10^^'' 
a.u. [4-9], the radiative and relativistic corrections have not been presented in the literature with required accuracy. 
Only recently high-precision variational calculations for the ro-vibrational states in the range of the total orbital 
momentum L = — 4 and vibrational quantum number u = — 4 for the and HD^ molecular ions along with 
relativistic and radiative corrections of orders RacO^ , RooOp' {ra j M) , RoaO^ , RaoO?{'mlM), and, partially, RooOt^ have 
been obtained [10]. 

In this paper we use the CODATA02 recommended values [11] for numerical calculations, and the atomic units 
{Ti = e = rrie = 1) are adopted. 

I. RADIATIVE CORRECTIONS OF ORDER R^a'^. 

For a given relative accuracy of ^^10^^*^ — 10"^^ recoil corrections of orders RooCt'^{iTi/M) and higher are small and 
may be neglected. That allows to reduce calculation of higher order corrections for the Coulomb three-body system 
to the problem of a bound electron in an external field. 

The radiative corrections of order RooO^ in the external field approximation are known in an analytic form [12, 13]: 
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The last equation includes both Dirac form factor and polarization operator contributions. 



II. RELATIVISTIC CORRECTIONS OF ORDER Ra^a^. 

The most problematic contribution of RoaO^ order is the relativistic correction for a Dirac electron. It can be obtain 
within the adiabatic two-center approximation as follows (for details, see [14]). 
We start from the nonrelativistic Schrodinger equation with the Hamiltonian: 

Ho = ^ + V, V = (2) 
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FIG. 1: Adiabatic "effective" potentials for the relativistic ma® order correction for molecular ion {Zi = Z2 = 1). Energies 
are in (atomic units) x a*. 




The total contribution to the energy of a bound electron at the RooCt'^ mgC^a^ order is defined by 
Ai;(6) = {HbQ{Eo - HoY^QHb) + (h^^^^ 



(3) 

Here iJ^^^ is the effective Hamiltonian for the interaction of an electron with the external field of two centers in this 
order, which can be expressed in a form: 
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Hb is the Breit-Pauli interaction: 
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Both terms in Eq. (3) are divergent. In order to remove the infinities a transformation to the second order term 
can be applied which separates a divergent part: 

(H's=Hb + {Ho -Eo)U + U{Ho - Eq) 

{HbQ{Eo - Ho)-^QHb) = {H'sQiEo - Ho)-^QH's) 

+{UHb+HbU) - 2 {U) {Hb) + {U{Ho - Eo)U) 

with U = ^[Z^/n + Z2M = -^V. 

The last three terms of the second equation in (6) can be recast in a form of a new effective interaction: 
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Here pi = ZiS{ri) and £i = —ZiVi/r^. 

Taking into account that ^'o is a solution of the Schrodinger equation Hq^q = Eq^q, one may obtain from the 
above the following finite expression [14]: 



A£(6) = {H'sQ{Eo - Ho)-'QH's) + + = {H'j,Q{Eo - H^r^QH'^) 

3Eo (V^) 5J?g {V) SwEq ((pi + p2)) (pVp) , {V){Hb) 
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This new expression along with modified second order iteration can be now calculated numerically. 

"Effective" potentials of IS.E^^\R) have been obtained for different bond lengths in [14]. Results are shown in 
Figure 1. 

Averaging them over the radial wave function of particular state one may get corresponding contribution to the 

energy of that state of order RooO^ ■ Results of numerical calculation of the relativistic corrections at this order are 
presented in Tables I and II. For the transition frequency this adiabatic approach provides about 3 significant digits. 
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TABLE I: Relativistic corrections of Roofna^ order (in units c'* x (1 a.u.)), Hj. 



V = V = 1 V = 2 11 = 3 V = 4: 



L = -0.042097 

L = l -0.042100 

L = 2 -0.042107 

L = 3 -0.042117 

L = 4 -0.042133 



0.042908 -0.043786 
0.042912 -0.043792 
0.042922 -0.043805 

0.042938 -0.043825 
0.042959 -0.043854 



0.044732 -0.045729 
0.044740 -0.045738 
0.044757 -0.045756 

0.044782 -0.045783 
0.044818 -0.045820 



TABLE IL Relativistic corrections of Rooma^ order (in units x (1 a.u.)), HD"*". 
V = t; = l 11 = 2 V = 3 V = 4 



L = -0.042043 

L=l -0.042045 

L = 2 -0.042050 

L = 3 -0.042058 

L = 4 -0.042069 



0.042738 -0.043483 
0.042741 -0.043487 
0.042748 -0.043496 
0.042759 -0.043510 
0.042773 -0.043529 



0.044278 -0.045126 
0.044284 -0.045132 
0.044295 -0.045146 
0.044312 -0.045167 
0.044336 -0.045195 



III. HIGHER ORDER RADIATIVE CORRECTIONS. 



The electron ground state wave function to a good extent may be approximated by 4>e{i^e) = C[V'is(ri) + '4>is{^2)], 
where tpis is the hydrogen ground state wave function. So, the most important RooO^ order contributions can be 
evaluated using this approximate wave function and the expressions: 
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where the constants Aqi, Aqq, and S50 are taken equal to the constants of the Is state of the hydrogen atom 
Aei = 5.419 . . . [15], Aqq = -30.924 . . . [16], and S50 = -21.556 . . . [17]. Worthy to say that the leading contribution 

{Rooa^ In^ a) is exact. 



IV. RESULTS AND CONCLUSION 



Various contributions to the frequency interval of the fundamental transition arc summarized in Table III. Uncer- 
tainty in orders RooCt^ and RooCx'^ arc primarily due to numerical uncertainty in calculation of leading order terms 
like (p^) in the Breit-Pauli Hamiltonian, or the Bethe logarithm, P{L,v) (see Refs. [18, 19] for the details), and can 
be improved by more extensive calculations. We estimate uncertainty due to finite size of nuclei as ~3 • 10~^ MHz for 



TABLE III: Summary of contributions to the {v = 0, L = 0) ^ (v' = 1, L' = 0) transition frequency (in MHz). 







HD+ 


AEnr 


65 687511.0686 


57349 439.9717 


AE^2 


1091.041(03) 


958.152(03) 


AE^3 


-276.544(02) 


-242.125(02) 


AE^A 


-1.997 


-1.748 


AE„5 


0.120(23) 


0.105(19) 


AEtot 


65 688 323.688(25) 


57350154.355(21) 
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these transitions. So, the latter is so far neghgiblc for the ro-vibrational spectroscopy. For the contribution of order 
Rooa^ the error bars are determined by the total contribution of the terms with coefficients A and B in Eq. (9). 

Recently, the {v, L) : (0, 2) (4, 3) ro-vibrational transition for the HD^ ion has been precisely measured in the 
experiment at the Dusseldorf university [20]. Comparison with theoretical calculation demonstrates a very good 
agreement: 

£;exp = 214 978 560.6(5) MHz 
Et^, = 214978 560.88(7) MHz 

In conclusion, the relativistic corrections of order RooOt^ allow to reduce the relative accuracy of the fundamental 
transition frequency in Hj to about 3 • 10~^ or 0.3 ppb. Further improvement we expect to achieve by numerical 
estimate of coefficients Aqi, Ago, and B50 from Eq. (9) using the two-center adiabatic (or external field) approximation. 
That may reduce the final uncertainty by a factor of 5-10 and the relative uncertainty to less than 10~^°. Eventually, 
it will make real the main goal of our studies: improving of the mp/me mass ratio from the ro-vibrational spectroscopy 
of H^ and HD+. 
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